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We point out inconsistencies in the recent paper by Oughstun et al. on Sommerfeld and Brillouin 
precursors [J. Opt. Soc. Am. B 27, 1664-1670 (2010)]. Their study is essentially numerical and, 
for the parameters used in their simulations, the difference between the two limits considered is not 
as clear-cut as they state. The steep rise of the Brillouin precursor obtained in the singular limit 
and analyzed as a distinguishing feature of this limit simply results from an unsuitable time scale. 
In fact, the rise of the precursor is progressive and is perfectly described by a Airy function. In 
the weak dispersion limit, the equivalence relation, established at great length in Section 3 of the 
paper, appears as an immediate result in the retarded-time picture. Last but not least, we show 
that, contrary to the authors claim, the precursors are catastrophically affected by the rise-time of 
the incident optical field, even when the latter is considerably faster than the medium relaxation 
time. 

OCIS codes: 260.2030, 320.5550, 320.2250. 



PACS numbers: 42.25.Bs, 42.50.Md, 41.20.Jb 

In a recent paper [lj, Oughstun et al. revisit the clas- 
sical problem of the propagation of a step modulated 
pulse in a Lorentz model medium. They specifically con- 
sider the case where the absorption line is narrow (sin- 
gular limit) and the one where the refractive index of 
the medium keeps very close to unity at every frequency 
(weak dispersion limit). The medium is characterized by 
its complex refractive index 
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Here to, u) p , ujq and S respectively designate the current 
optical frequency, the plasma frequency, the resonance 
frequency and the damping or relaxation rate. The wave 
propagates in the z-direction. In the following we use the 
retarded time t, equal to the real time minus z/c where 
c is the velocity of light in vacuum (retarded time pic- 
ture). The medium is then characterized by the transfer 
function 



H{uj) = cxp i — z (n(uj) 
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and the field transmitted E{z, t) at the abscissa z reads 



E(z,t) = — 
Ztt 
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H(oj)E(0,cjj)exp(-iut)duj. (3) 
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Here a is a positive constant and E(0,uj) is the Fourier 
transform of the incident field E(0,t). In the latter 
is assumed to have the idealized form 



E(0,t) = Q(t) sin(u c t) 



(4) 



where 0(i) is the unit step function and ui c is the fre- 
quency of the optical carrier. 



Although it abundantly refers to the theoretical results 
obtained by the asymptotic method, the study reported 
in [l| is mainly numerical. All the simulations are made 
for lu q = 3.9 x 10 14 rad/s [2] and uj c = 3.0 x 10 14 rad/s in 
a normal dispersion region. The singular and weak dis- 
persion limits are respectively attained when 5 <iC loq and 
ujp <C Sujo [see Eq. ((TJ]. Oughstun et al. emphasize that 
these two limiting cases "are fundamentally different in 
their effects upon propagation" but, surprisingly enough, 
they take for their simulations in the weak dispersion 
limit a value of 6 for which the singular limit nearly 
holds (5 < wo/100). Consequently, mutatis mutandis, 
the results obtained in the two limits appears qualita- 
tively similar. The steep rise of the Brillouin precursor 
obtained in the singular limit (their Fig. 3) and analyzed 
as a distinguishing feature of this limit is only due to an 
unsuitable time scale. As shown in our FigJTJ obtained 
for the same values of the parameters, the rise is quite 
progressive and well reproduced by a Airy function. This 
result, established by Brillouin himself in 1932 [3j, is eas- 
ily retrieved from Eqs (JTJ) , (J2J) , and ([3]) . Anticipating that 
the beginning of the precursor involves frequencies u) such 
that 8 <C u <C u) C) we use the approximate relations 
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Introducing the new retarded time 
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and E(0,w) « 1/cj c 
t' = t — % with % = 
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function then reads as exp [iw 3 / (36 3 )] where 
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Figure 1: Brillouin precursor obtained for the parameters of 
the figure 3 of PJ, that is for ojq = 3.9 x 10 14 rad/s, iv c — 
3.0 x 10 14 rad/s, ui p = 3.05 x 10 14 rad/s , 5 = 3.02 x 10 10 rad/s, 
and z — 7.232 x 10 _2 m. The full line is the exact numerical 
solution obtained by fast Fourier transform (FFT) and the 
dashed line is the approximate analytical solution given by 
Eq[7l 



We finally get 



E(z,t') 
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exp ( — T - iut' ] dw = — Ai(-bt') 
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where Ai(x) is the Airy function [4|. As it appears Fig[TJ 
this analytic expression perfectly fits not only the rise of 
the precursor but also a significant number of its subse- 
quent oscillations. Equations ([6]) and ([7]) also provide a 
simple evidence of the z -1 / 3 dependence of the precursor 
amplitude on the propagation distance. Figure [2] shows 
the Brillouin precursor obtained in the conditions of the 
Figure 5 of intended to illustrate the specific case of 
the weak dispersion limit. As indicated previously, the 
conditions of the singular limit are approximately met 
and this explains why the first oscillation of the precur- 
sor and thus its peak amplitude are well reproduced by 
Eq. (J7|). Note that the rise-time of the precursor, pro- 
portional to 1/&, is 6 times faster than in the previous 
case. 

In their study of the weak dispersion limit, Oughstun 
et al. [l[ mention a "curious difficulty in the numerical 
FFT simulation of pulse propagation" and, in order to 
overcome it, they develop at great length an equivalence 
relation. In fact the difficulty is completely avoided in the 
retarded-time picture and their equivalence relation then 
appears as an immediate result. In the weak dispersion 
limit, the transfer function, as given by Eq.©, is indeed 
reduced to 



H{uj) exp 




(8) 




This only expression shows that, for given ujq and 5, all 
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Figure 2: Same as Fig.l for the parameters of the figure 5 of 
[U, that is for w = 3.9 x 10 14 rad/s, cj c = 3.0 x 10 14 rad/s, 
uj p = 3.05 x 10 12 rad/s , 6 = 3.02 x 10 12 rad/s, and z = 2.290m. 
The rapid oscillations of small amplitude superposed upon 
the beginning of the Brillouin precursor are the end of the 
Sommerfeld precursor. 



the media having the same uj^z (and thus the same opti- 
cal thickness at any reference frequency) are equivalent. 

In a real experiment, the incident field is obviously 
turned on in a finite time. Oughstun et al. state that 
their "results will remain valid for a non instantaneous 
turn-on signal provided that the signal turn-on time T r is 
faster than the characteristic relaxation time 1/(5". This 
is grossly false. Figure [3] shows the results of a simula- 
tion performed for the parameters of their figure 3 and 
10 - 90% turn-on times (a) T r = 0, (b) T r = r/500, 
and (c) T r — r/100 where r = 1/6 is the relaxation or 
damping time (r « 33ps) . The rise of the incident field is 
modelled by replacing the step- function Q(t) appearing 
in the idealized form of Eq. (0} by f(t) = \ [1 + erf(£i)] 
where erf (a;) designates the error function. The corre- 
sponding rise-time is T r ss 1.8/£ and f(t) — > <d(t) when 
£ — > oo. We see that the effect of the rise-time is actually 
catastrophic. For a rise-time as fast as r/500 (« 66fs 
), the Sommerfeld precursor is absent. The maximum 
of the Brillouin precursor is significantly time-delayed 
and its amplitude Ag is reduced by a factor of about 
300, becoming comparable to the amplitude e -10 of the 
"main field". The reduction factor obviously depends on 
the optical thickness and, consequently, the power law 
Ab oc z -1 / 3 , obtained for T r = 0, breaks down. Finally 
the curve (c) of FigJ3] shows that both precursors prac- 
tically vanish for T r — r/100 . By means of other sim- 
ulations, we find that the Sommerfeld precursor is very 
attenuated as soon as T r > r/5000 « 2/u c and that a 
correct reproduction of both precursors as obtained for 
T r = requires that T r < 1/lo c . This condition results 
from the fact that the Sommerfeld (Brillouin) precursor 
mainly involves frequencies uj 3> uj c (lo <C uj c ) and thus 
is excited by the corresponding frequencies contained in 
the spectrum of the incident field. From the expression of 
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Figure 3: Effect of the rise-time of the incident field on the 
transmitted field. The parameters are those of the figure 3 
of [H, that is w = 3.9 x f0 14 rad/s, uj c = 3.0 x f0 14 rad/s, 
L0 P = 3.05 x f0 14 rad/s , S = 3.02 x f0 10 rad/s, and z = 7.232 x 
10 _2 m. The vertical dashed line indicates the group delay 
T s = \u=w c where tp is the argument of H(u>) (r 9 = 471ps). 
T g fixes the arrival of the "main field" whose amplitude, equal 
to exp (— 10) ~ 4.5 x 10~ 5 , obviously does not depend on the 
rise-time. The different curves are obtained for (a) T r = 0, 
(b) T r = r/500, and (c) T r = r/100 where r = 1/8 is the 
relaxation or damping time of the medium. 



f(t), it is easily shown that the effect of a finite rise-time 
T r (T r w 1.8/C ) is to divide E(Q,u) by exp [w 2 / (4£ 2 )] 



when ui/uj c —> oo and by exp [w 2 / (4£ 2 )] when lu/lj c —5- 0. 
These expressions explain why the Sommerfeld precursor 
is much more affected by the rise-time effects than the 
Brillouin precursor and enable one to predict that, for 
T r = 1/cJc (£ ~ l-8w c ), the amplitude Ab of the Bril- 
louin precursor will be about 7.5% below that obtained 
for T r — (result confirmed by an exact numerical cal- 
culation). 

The observation of a signal close to that shown on 
the figure 3 of [l| requires not only that the rise-time 
of the incident field does not exceed 1/oj c (3.3fs) but also 
that its subsequent amplitude remains nearly constant 
during a time that would be more than four orders of 
magnitude longer. The fulfilment of this double condi- 
tion, either with a pulsed laser or with a continuous wave 
laser followed by a modulator, appears to be quite unre- 
alistic from an experimental viewpoint. More generally, 
due to such temporal constraints, we don't see how the 
Brillouin precursor could be actually used in optics as 
a tool for imaging through a dense medium, opaque at 
the carrier frequency. Anyway, the problem of optical 
precursors is of fundamental interest from a theoretical 
viewpoint, even if the experiment imagined by Sommer- 
feld and Brillouin may appear as a gedankenexperiment. 
In this spirit, we are examining what the precursors be- 
come when the incident field E(0, t) = 9(t) sin(u; c £), con- 
stantly considered in the theoretical papers, is replaced 
by E(0, t) = 0(i) cos(w c i). The true discontinuity of the 
incident field at the initial time then leads to radically 
new effects that we will discuss in a forthcoming paper. 
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